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ABSTRACT 

In this paper we derive an improved core mass function (CMF) for the Pipe 
Nebula from a detailed comparison between measurements of visual extinction 
and molecular-line emission. We have compiled a refined sample of 201 dense 
cores toward the Pipe Nebula using a 2-dimensional threshold identification al- 
gorithm informed by recent simulations of dense core populations. Measurements 
of radial velocities using complimentary C^^O (1-0) observations enable us to cull 
out from this sample those 43 extinction peaks that are either not associated with 
dense gas or are not physically associated with the Pipe Nebula. Moreover, we 
use the derived C^^O central velocities to differentiate between single cores with 
internal structure and blends of two or more physically distinct cores, superposed 
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along the same line-of-sight. We then are able to produce a more robust dense 
core sample for future follow-up studies and a more reliable CMF than was possi- 
ble previously. We confirm earlier indications that the CMF for the Pipe Nebula 
departs from a single power- law like form with a break or knee at M ~ 2.7 ± 
1.3 Mq. Moreover, we also confirm that the CMF exhibits a similar shape to the 
stellar IMF, but is scaled to higher masses by a factor of ~ 4.5. We interpret 
this difference in scaling to be a measure of the star formation efficiency (22 ± 
8%). This supports earlier suggestions that the stellar IMF may originate more 
or less directly from the CMF. 

Subject headings: stars: formation-dust, extinction-ISM: globules-ISM: molecules- 
stars: luminosity function, mass function 



1. Introduction 

Recent studies of dense cores suggest that the fundamental mass distribution of stars 
may be set during the early stages of core formation. This is primarily because of the 
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If the stellar IMF is in fact predetermined by the form of the CMF, then the origin of the 
stellar IMF may be directly linked to the origin of dense cores. Thus, understanding the 
connection between the CMF and the stellar IMF is essential for theories of star formation. 

Power-laws are by nature scale-free and the apparent similarity of the CMF and IMF 
slopes doe s not necessarily imply a one-to-one correspondence between the CMF and the 
IMF (e.g.. lSwift fc Williams! I2OO8I ). However, because the shape of t he stellar IMF flattens 
and departs from a single power-law at low masses (e.g.. lKroupall200ll ). any similar flattening 
measured for a CMF would strengthen the idea of a very direct connection between the two 
distributions. Moreover, the characteristic mass, that is, the mass at which the distribution 
departs from a single power-law form, provides a deflnite physical scaling for a mass function 
and differences in the characteristic masses of the CMF and IMF would be a measure of the 
efficiency of star formation (e.g. 



Alves et al.l 120071 ) 



The reliability of any detailed statements made about the shape of the CMF, or its 
connection to the stellar IMF, depends sensitively on the uncertainties involved in the gen- 
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eration of the CMF. Small number statistics, core crowding and the (in)accuracy of the 
core masses can very easily iiitrodu ce significant uncertainties in the derived shape of the 



CMF (e.g., iKainulainen et al.ll2008l ). Even the operational definit ion of a core caii be a 



troubling source o f uncertainty in the determination of a CMF (e.g., ISwift &: Williams! 12008 



Smith et all 120081 ) . 



To minimize such uncertainties and to obtain a meaningful measurement of a CMF 
we first require a large sample of dense cores with accurate measurements of their proper- 
ties. Ideally we desire a sample of starless cores in a single molecular cloud. Fortunately, 
the Pipe Nebula, at a distance of 130 pc, is well suited for this purpose. The cloud ex- 
hibits little evidence for active star formation and a detailed dust extinction map of the 
comple x exists containing nunierous core-like st ructures w hose masses c an be precisely mea- 
sured dhombardi eTaP bood : klves et al.lboOTh . Indeed, klves et al.l J2007h identified 159 
cores with masses between 0.2-20.0 Mq in the Pipe Nebula. Subsequent molecular-line ob- 
servations demonstrat ed that many of these were in fact dense cores (n(H2) > 10^ cm~'^; 



Rathborne et al.ll2008l ) . The CMF derived for this core population was not a single power- 
law function. Instead, the CMF exhibited a clear break at a mass of ~ 2 M©, suggesting 
that an overall efficiency of about 30% would characterize the star formation process if these 
cores evolve to make stars (lAlves et al.l 120071 ). 



Kainulainen et al.l (120081 ) performed simulations of the Pipe core population and inves- 



tigated the sensitivity of the deriv e d CM F to the details of the core extraction algorithm 
that was employed by lAlves et al.l (120071 ). The simulations confirmed the presence of the 
break in the derived CMF and revealed that the mass range where the break occurs is not 
affected by either incompleteness or biases introduced by the extraction algorithm. However, 
they found that the fidelity of the overall shape of the CMF did depend on the choice of 
input parameters of the core extraction algorithm and they determined optimu m values for 



the se arch parameters. These parameters differed somewhat from those used by lAlves et al. 



(120071 ). Moreover, the simulations also indicated that extractions sometimes produce spuri- 
ous cores. Finally there are also uncertainties inherent in the use of 2-dimensional data such 
as extinction maps, to accurately define a core population. In particular, such data cannot 
distinguish between any foreground or background extinction features or cores. In addition 
it is difficult to distinguish physically related substructure within a single object from a blend 
of physically distinct cores along a similar line-of-sight. These considerations suggested that 
a re-analysis of the Pipe CMF was warranted, especially given the ramifications of the results 
for core and star formation. 



To derive a more reliable CMF for the Pipe Nebula and investigate its relation to the 
stellar IMF, we have undertaken a study of cores within the Pipe Nebula using a combination 
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of visual extinction and new C^^O (1-0) molecular line data. Because the core properties 
are derived directly from how cores are defined and extracted, one needs to pay careful 
attention to how this is done. Here we use the visual ext i nction image and a 2-dimensional 
clump finding algorithm, guided by the lKainulainen et al.l (120081 ) simulations, to identify and 
extract discrete extinction features as candidate cores. The C^^O (1-0) emission is then used 
to determine if the candidate cores are (1) associated with dense gas, (2) associated with 
the Pipe Nebula, and (3) separate or physically related structures. Because 2-dimensional 
finding algorithms can pick up spurious or unrelated cores, we use C^^O (1-0) emission to 
filter the sample to include only real dense cores that are associated with the molecular 
cloud. This is particularly important because the inclusion of spurious or unrelated cores, as 
well as the artificial merging or separation of extinction features, will significantly affect the 
shape of the CMF and any conclusions derived from it. With a reliable CMF we can then 
make detailed comparisons between it and the stellar IMF to gain a better understanding of 
the connection between dense cores and the star formation process. 



2. Identification of discrete extinction peaks 



The identification of discrete extincti on peaks within t he Pip e Nebula was performed 
using the visual extinction [Ay) image of iLombardi et al.l (120061 ). This image was gener- 
ated using the NICER method which utilized the JHK photometry of over 4.5 million stars 
within the 2MASS catalog. In addition to the many compact extinction features, this im- 
age also contains significant and variable extinction from the lower-column density material 
associated with the molecular cloud. To reliably extract discrete extinction peaks a back- 
ground subtraction was performed using a wav elet decomposition which filters out structures 
larger than the specified size scale (0.30pc; see lAlves et al.ll2007l for details). This procedure 
produces a smoothly varying background extinction image and a 'cores-only' image. The 
extinction features were extracted from the cores-only, background-subtracted extinction im- 
age. All core parameters (coordinates, peak Ay, radius, and mass) were measured directly 
from this cores-only map. 

To identify the discrete extinction peaks, we use the 2- dimensional version of the auto- 
mated algorithm clumpfind (clfind2d; IWilliams et al.lll994j ). Clfind2d searches through an 
image using iso-brightness surfaces to identify contiguous emission features without assum- 
ing an a priori shape. The iso-brightness surfaces are identified through a series of contour 
levels. In contrast to the 3-dimensional version of the algorithm, clfind2d allows users to in- 
put arbitrary levels to determine the contouring. The algorithm starts its search for emission 
features at the pixel with the peak brightness in the image and steps down from this using 
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the user-defined contour levels. During each iteration, the algorithm finds all contiguous pix- 
els between each particular contour level and the next level down. If the contiguous pixels 
are isolated from any previously identified emission features, then they are assigned to a new 
clump. If they are connected they are assigned to a pre-existing clump. For cases of isolated 
emission features, this procedure is straightf orward. However, in t he case of blended emission 
features, a 'friends-of-friends' algorithm (see I Williams et al.lll994j for more details) is used to 
separate the emission. The algorithm iterates until the lowest contour level is reached. At 
all levels clfind2d requires that each new clump is greater than a specified number of pixels. 
To be conservative, we set this minimum number of pixels to be 20. Assuming this area is a 
circle, then this size corresponds to 2.5 times the angular resolution of the visual extinction 
image. 



The original list of cores toward the Pipe Nebula was compiled by lAlves et al.l (120071 ) us- 
ing the same wavelet decomposition and clfind2d algorithm described above. They identified 
159 cores using extinction contour levels of 1.2, 4.0, and 6.0 mags. Because the extinction 
contouring was truncated at 6 magnitudes, many of the higher extinction regions were con- 
sidered to be single extinction features when in fact some consist of multiple well separated 
extinction peaks. 

To improve on the previous list of cores and to identify them in a homogeneous way at all 
emission levels, we choose to define the contour levels input into clfind2d using discreet values 
of the small scale variation in the background [ctav)- This was estimated by calculating the 
standard deviation in the extinction values in ~1 pc^ regions across the background map. 
Because the extinction which results from the larger-scale molecular cloud varies across the 
Pipe Nebula, we use the value of aAy to be 0.4 mags . This value is higher than the individual 
noise in each pixel (0.2 mags; iLombardi et al.ll2006l ). 



Our choice of clumpfind input parameters is based on the recent simulations of lKainulainen et al, 

(120081 ). To best match the real data, these simulations placed elliptical Gaussian cores of 



various masses on the background image of th e Pipe Nebula using the observed median core 
separation and column density distribution of lAlves et al.l (120071 ). Similar to the real data, 
a background subtraction was performed using a wavelet decomposition. The extraction of 
cores was then performed on the background-subtracted image using a range of clumpfind 
contour levels. These simulations were investigating the effect different clumpfind contour 
levels have on the der ived core properties and the recoverability and completeness of the 
resulting CMF. While iKainulainen et al.l (120081 ) find that the exact choice of contour levels 
has little effect on the overall mass completeness limit, the derived CMF is slightly more 
complete when using clumpfind contours that include the extinction values down to the im- 
age noise. More importantly, however, these simulations show that the degree of crowding 
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within a molecular cloud can significantly effect both the measured core parameters and the 
derived CMF. 

The simulations reveal that the Pipe CMF is 90% complete to a mass of 0.5 Mq when 
using contour levels starting at 0.4 mag {laAv) increasing in steps of 1.2 mag (3cr^^). 
In a typical simulation where ~230 cores were input, these clumpfind parameters extracted 
~190 cores. While the number of extracted cores is close to the number input, a visual 
inspection of the positions of the input and extracted cores reveal significant differences. We 
found that ~60% of the extracted cores match well with those input, while ~33% of the 
input cores were not extracted. More concerning, however, is the large number of spurious 
cores that were extracted: ~ 40% of those extracted did not match an input core. Because 
these are typically low mass cores, it is likely that the majority of these spurious cores may 
have been introduced by the wavelet decomposition and are not real dense cores. Moreover, 
the simulations also reveal that the process of wavelet decomposition filtered out some of 
the input cores. While these effects are most obvious for cores that are below the mass 
completeness limit, the use of such a low contour level artificially merges some of these noise 
fiuctuations with 'real' cores, thereby increasing their mass which will then alter the shape 
of the resulting CMF. 

To reduce the number of these spurious cores extracted from the real data, we chose 
to start the contouring at 1.2 mag (3(1^^,) rather than 0.4 mag {Ictav)- Specifically, we use 
contour levels starting at 1.2 mag {?>aAy) and increasing in steps of 1.2 mag (3o"a^), i.e. 1.2, 
2.4, 3.6, 4.8, 6.0, 7.2, 8.4, 9.6, 10.8, 12.0, 13.2, and 14.4 mags. 



With these contour levels iKainulainen et al.l (120081 ) find that from twenty realizations of 



the simulation, the mean number of cores extracted is ~174 (when ~240 are input). From 
the extracted cores in a single simulation, ~ 150 of them are coincident with an input core, 
while ~ 26 of them have no corresponding input core and are, thus, spurious detections. The 
simulations reveal that the mass completeness limit is ~0.65 Mq. Most (97%) of the input 
cores that were not extracted were below this mass completeness limit. Moreover, the most 



massive core missed was typically not greater than ~ 1.5 Mq. In addition. IKainulainen et al. 



( 120081 ) find that only about 15% of the input cores are blended with their nearest neighbor. 



Using contour levels starting at 1.2 mag (3cr^^) and increasing in steps of 1.2 mag (3(7^^) 
on the observed visual extinction cores-only image, we extract 201 discrete extinction peaks 
from the real data0. We use this list for all further analysis. 



^These parameters appear to characterize the underlying data weU, the exception being toward the 
densest region in the Pipe Nebula, Barnard 59. Using these parameters clfind2d identifies 6 extinction peaks 
toward Barnard 59. Because of the high extinction and incomplete number of stars toward this region in the 
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3. Molecular line observations 

Observations of C^'^O (1-0) molecular line emission toward all 201 extinction peaks 
identified toward the Pipe Nebula were obtained using the 12 m Arizona Radio Observatory 
(ARO) and the 22 m Mopra telescopes. The ARO C^^O (1-0) observations were conducted 
over two periods in 2005-2006. In total, 102 of the 201 extinction peaks were observed with 
the ARO. Ninety -four of these positions corresponded to cores wit hin the origina l list o f 



Alves et al.l (120071 ). The details of these observations are presented in lMuench et al.l ( 120071 ). 



Spectra were obtained with the Mopra telescope toward all remaining identified extinc- 
tion peaks for which we had no previous C^'^0 (1-0) data. In total, an additional 99 positions 
were observed over the periods 2007 July 3-8 and 2008 July 20-25. In all cases the spectra 
were obtained toward the peak extinction identified using clfind2d and the 2MASS visual 
extinction image {Qhpbw ~ 1')- 

For the Mopra observations the 8 GHz spectrometer MOPS was used to simultaneously 
observe the transitions ^^CO, ^^CO, C^'^O, and C^''0 (1-0). The spectrometer was used 
in 'zoom' mode such that one spectral window covered each of the lines of interest. Each 
window was 137.5 MHz wide and contained 4096 channels in both orthogonal polarizations. 
This produced a velocity resolution of 0.09 kms~^ for t he C^'^0 (1-0) transition. At these 



frequencies the Mopra beam is ~ 33" (ILadd et al.ll2005l ). 



All spectra were obtained as four 5 minute integrations in the position switched mode. A 
common 'off' position was used in all cases and was selected to lie within a region of the Pipe 
Nebula which is devoid of cores and molecular gas (a!=17:22:51.75, 5=— 25:24:12.39, J2000). 
The telescope pointing was checked approximately every hour using a suitably bright, nearby 
maser. 

The spectra were reduced using the ATNF Spectral Analysis Package (ASAP) and were 
initially baselined subtracted before averaged using a system temperature (Tsys) weighting. 
The typical Tgys for these observations was ~ 320 K. Gaussian profiles were fit to each 
spectrum to determine its peak temperature (TX), central velocity (vlsr), line-width (AV), 
and integrated intensity (1). To be considered a significant detection, we require that the 
integrated intensity be greater than 3 times the rms noise. All raw data are in the 
scale. The final spectra have a typical rms noise of ~0.02 K channel"^. To convert to 



2MASS image, the extinction image has artificial structm'e which clfind2d breaks up in to muhiplc peaks. 
Higher-angular resolution images show, however, that this core though mostly smooth may in fact contain 
substructure consisting of 2-3 extinction peaks at its center. For this work, however, we incorporate all the 
extinction within this region in to a single core. 
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th e main beam brigh tness, Tmb, one needs to use the main beam efficiency of 0.43 as hsted 



m 



Ladd et ahl (120051 ) . Although the spectra have a velocity resolution of ~0.09 kms ^, the 
typical error on the measurement of vlsr from the Gaussian fits were <0.01 kms~^. 



4. Results 



Having identified a large number of extinction peaks within the background-subtracted 
extinction image, we now use the C^^O (1-0) emission to determine if they are associated 
with dense (~10^cm~^) gas and which are associated with the Pipe molecular cloud. Exclud- 
ing any spurious extinction peaks that may have been artificially included via the wavelet 
decomposition and any that are unrelated to the region is crucial in generating a reliable 
CMF. In this section we also use the C^^O (1-0) emission to help determine whether nearby 
extinction features are physically separate or part of the same physical core. With these dis- 
tinctions we can then derive reliable core ensemble properties such as radius, mass, density, 
and non-thermal line-width. 



4.1. Dense gas associated with the extinction peaks 

Although we also have ^^CO and ^^CO toward a large number of the extinction peaks 
within the Pipe Nebula, we choose to focus on the C^^O emission. The critical density for 
excitation of C^^O (1-0) is 6xl0^cm~^ and because the isotope is rarer and the emission 
is optically thin, it is a far superior tracer of dense gas within a core than either ^^CO or 
^^CO. Combining the ARO and Mopra C^^O (1-0) observations, we find that 93% of the 
extinction peaks have significant C^^O molecular line emission. Thus, we find that 188 of 
the extinction peaks are associated with dense gas. 

The extinction peaks with no detected C^^O (1-0) emission have low peak extinctions 
{Ay < 3.6 mags) and low masses (M < 1.3 Mq). Moreover, the majority of these are isolated 
features and are located toward the edges of the extinction image. It is likely that these 
extinction peaks are not real dense cores and could simply be noise fluctuations artifically 
added in by the wavelet decomposition. 



4.2. Association with the Pipe Molecular Cloud 

Previous large-scale ^^CO observations show that the main body of the Pipe Nebula is 
associated with molecular line emission in the velocity range of 2 kms~^< vlsr< 8 kms~^ 
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( lOnishi et al.lll999l ). Figure [T] shows the measured C^^O central velocities [vlsr) toward all 
the extinction peaks with significant C^^O emission. We find that the emission from the 188 
extinction peaks falls within a narrow range. A Gaussian fit to the distribution reveals that 
it is centered at a vlsr of ~ 3.9 kms~^ with a FWHM of ~ 2 kms~^. 

To determine which extinction peaks are associated with the Pipe Nebula, we use the 
range of 1.3 kms~^ < vlsr < 6.4 kms~^ (center vlsr =t 3 times the FWHM). Of the 188 
extinction peaks that are associated with C^^O emission, we find that 158 of them have 
velocities in this range. Thus, we will consider only these extinction peaks as associated 
with the Pipe molecular cloud. 

We find that 30 of the extinction peaks have associated molecular gas that differs sig- 
nificantly from this velocity range, however. This gas, at velocities of ~— 5 to kms~^ and 
~15 to 20 kms~^, is likely associated with foreground and/or background molecular clouds. 
Indeed, there are 16 extinction peaks with a vlsr of ~ 20 kms~^, that are localized in 
Galactic longitude and latitude {i ~3°, 6 ~ 4.2°). These extinction peaks are most likely 
related to a common background molecular cloud. 

For all further analysis, we consider only the 158 extinction peaks that have dense gas 
at velocities associated with the Pipe molecular cloud. 



4.3. Multiple, nearby extinction peaks: superimposed cores or cores with 

substructure? 

If cores are centrally condensed objects that will give rise to a single star, then their 
extinction profiles should refiect their density gradients and be approximately symmetric, 
centrally peaked, and well defined above the background. While many isolated extinction 
features within the Pipe Nebula have these simple characteristics and can be identified as 
discrete cores, there are several extinction features that show more complex structure. 

Because of the potential infiuence of crowding, many of these complex extinction struc- 
tures may actually comprise multiple cores. With the 2-dimensional visual extinction image 
alone, it is impossible to determine if these complex structures comprise individual cores that 
are superimposed along the line of sight or if they are single cores with substructure that may 
not necessarily give rise to separate star-forming events. Because we are using 2-dimensional 
extinction data to determine the core properties, we need to distinguish between these sce- 
narios. An incorrect extraction of the core parameters will directly infiuence their derived 
masses and, hence, the shape of the resulting CMF. For example, the erroneous merging 
of many extinction features into a single structure will artificially increase the number of 
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high mass cores and decrease the number of low mass cores. In contrast, if the substructure 
within a larger core is artifically separated into many individual objects, then there will be 
a deficiency of high mass cores and an abundance of lower mass cores. In either case, the 
derived CMF will be significantly affected. 

Because we are interested in deriving reliable masses and, hence, the CMF for cores 
within the Pipe Nebula, we need to distinguish between extinction features that are the 
superposition of physically unrelated cores from those that may have internal substructure. 
As a first attempt at distinguishing between these scenarios, we use the vlsr derived from 
the C^^^O (1-0) emission. Because most (~80%) of the extinction peaks are well-defined 
isolated features, it is only necessary to consider those larger extinction features that appear 
to contain substructure. Thus, for this analysis, we consider only multiple extinction peaks 
that happen to lie within a larger structure. 

We define the larger extinction structures as contiguous regions of visual extinctions 
greater than 3cr (> 1.2 mags) in the background-subtracted image. Given the clumpfind 
contouring levels described in § [21 we find that 28 large extinction structures contain multiple 
extinction peaks. Because we have measured C^^O (1-0) emission toward all the extinction 
peaks, we can calculate the difference in the vlsr of each peak with respect to all other peaks 
within the larger structure B If the velocity difference (5V) between two extinction peaks 
is greater than the 1-dimensional projected sound speed in a 10 K gas (ci^ = 0.12 kms~^), 
then we assume that the C^^O (1-0) emission is arising from physically different cores that 
happen to be adjacent because they are superimposed along the line of sight. On the other 
hand, if the velocity difference is less than cid, then the extinction peaks may in fact be 
physically associated. To be considered part of the same physical core, we also require that 
the separation between the extinction peaks with velocities less than cid is smaller than its 
Jeans length, Rj. The Jeans length was calculated via the expression 

y S-KGrUpp 

where k is the Boltzmann constant, Tq is the gas temperature, G is the gravitational constant, 
rrip is the mass of a proton, and p is the mass density. In all cases we assume Tg=10K. 



^The molecular line emission associated with the isolated extinction peaks all correspond to a si ngle Gaus- 
sian C^^O (1-0) line profile. Some positions within the Pipe molecular ring (jMuench et al.l 120071 ). however, 
have multiple C^*0 (1-0) line profiles. In these cases we consider all emission lines when determining the 
relative vlsr between neighboring extinction peaks. 
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Figure [2] shows three examples of how we determine whether adjacent extinction peaks 
are physically associated. The left panels show the background-subtracted extinction images 
with contours as defined in § |5J The crosses mark the position of each extinction peak 
identified. The right panels show the resulting cores after taking into account the velocity 
differences and Jeans lengths of the extinction peaks. Marked on these images are the 
central velocities determined from the C^^O (1-0) emission (V), the distance to the nearest 
extinction peak (D), and the Jeans length (Rj). Using the criteria, 5V < Cic and D < Rj, 
we have determined if adjacent extinction peaks are isolated entities or part of the same 
physical core. In these images the color scale represents the area that is assigned to each 
core. In some cases we find that highly non-symmetric, complex extinction features have 
almost the same central velocity: these cores appear to have low-extinction tails (e.g. top 
and middle panels in Fig. [2]). In other cases, however, we find that adjacent extinction peaks 
within the same large scale extinction feature can have very different velocities (e.g. lower 
panels of Fig. [2]). We assume that these are physically differentiated or distinct cores. 

This procedure is summarized in Figure [3] which plots the velocity difference (^V) be- 
tween each extinction peak and every other within each large extinction structure against 
the distance to the extinction peak's nearest neighbor (D). The filled circles mark the 6V 
and D for extinction peaks that were determined to be part of the same core. In total, 
we find that 41 extinction peaks have 6V < cid and D < Rj. These were merged into 17 
cores for inclusion in the final list. We note, however, that this method will only allow us 
to separate cores that have significant relative motion along the line-of-sight. If the relative 
motion between adjacent cores is in fact in the plane of the sky, we will h ave incorrectly 



merged them. Considering a core to core velocity dispersion of ~ 0.4 kms ^ (IMuench et al. 



20071 ). we crudely estimate that no more than 12 of the 41 cores may have been incorrectly 



merged. 

To determine the masses of the merged cores, we calculate the sum of the total extinction 
within the area associated with each individual extinction peak. In all cases, we derive sizes 
using the total number of pixels, converting this to a radius assuming that the total area 
is a circle. For the remaining extinction peaks, we assume that they are unrelated cores. 
Although many of these extinction peaks have a neighbor that is either close in velocity or 
distance, they do not satisfy both of these criteria. Thus, we assume these are simply chance 
superpositions of unrelated cores along the line of sight. To determine their masses we use 
the total extinction directly output from clumpfind. 

Thus, after consideration of the velocity differences and Jeans lengths of the identified 
extinction peaks, we find that there are 134 physically distinct dense cores associated with 
the Pipe Nebula. Figure H] shows the location and approximate extent of the cores identified 
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within the Pipe Nebula. 



4.4. Core Mass Functions 

The derived CMFs are represented in Figure [5] a s binned histograms . Included on each 



plot is a scaled stellar IMF for the Trapezium cluster (IMuench et al.ll2002[ ) . The vertical dot 



ted lin e marks the mass completeness limit calculated in the simulations of iKainulainen et al. 



(l2008h 



The four panels in this figure represent the different stages in the selection process 
for dense cores outlined above and illustrate the changes in the shape of the CMF that 
occur as the result of the refinement of the core sample as one includes more information 
from the C^^O (1-0) emission. Figure [5] (a) shows the derived CMF for the 201 extinction 
peaks identified from the extinction image using clfind2d and the contour levels described 
in § 121 This represents the CMF derived using a blind application of clfind2d and includes 
all extinction peaks as cores, regardless of whether or not they have associated C^^O (1-0) 
emission and, thus, dense gas. This will include spurious and unrelated background cores 
and may also artifically break high-mass cores into many lower-mass objects. 

If we then consider only extinction peaks that have C^^O (1-0) emission and dense 
gas, the number of extinction peaks included in the CMF is reduced to 188, the shape of 
which changes only slightly for the lowest masses (M < 2 Mq; Fig. [5] (b)). This is not 
surprising considering the identification of spurious cores is a significant effect only at the 
lowest extinction levels. Similarly, our exclusion of cores lying outside the radial velocity 
range considered for the Pipe Nebula (30 cores; Fig. \5\ (c)) slightly alters the shape of the 
CMF for masses < 2 Mq. 

A more significant change in the shape of the CMF occurs when we consider the C^^O (1- 
0) velocity differences and angular distances between neighboring extinction peaks, as shown 
in Figure (d). In this case, extinction features are merged if their velocity differences are 
less than the 1-dimensional projected sound speed and separations are less than a Jeans 
length. This results in an increase in the number of cores at the high-mass end (M > 5 M©) 
and a decrease in the number of cores between masses of ~0.3 and ~3 Mq. 

Regardless of the exact shape, it appears that none of the CMFs shown in Figure 
are characterized by a single power-law. Instead, there is a break from a single power-law 
form providing a physical scale or characteristic mass for the CMF around ~ 2-3 Mq. The 
lowest mass bins in the CMF are most likely seriously effected by incompleteness. Thus, 
the position of the peak and the turnover at the very lowest masses (M~ 0.4 Mq) may be 
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unreliable. The simulations show that this mass range is most affected by the incompleteness 
due to the wavelet decomposition. While all the CMFs peak at roughly the same mass and 
show a break that is well above the mass completeness limit, the most likely to be reliably 
tracing the underlying distribution of core masses is the CMF shown in Figure O (d). Thus, 
from here forward we adopt this distribution as the CMF for the Pipe Nebula. 

It is of interest to core and star formation studies to compare the forms of the CMF and 
stellar IMF. To achieve this quantitatively, we have performed a minimization between 
the Pipe CMF and the stellar IMF by simultaneously scaling the IMF in both the x and 
y directions. The scaling factor in x direction will give the mass scaling between the CMF 
and the IMF. Assuming that each core will give rise to a ~1 star on average, this offset 
will give an estimate of the star- formation efficiency (SFE), that is, how much of the typical 
core mass is converted into the final stellar mass. We have also calculated the probability 
and use this to estimate the errors. The quoted errors are calculated from the range in the 
values for which the ^ probability is greater than 95 % (i.e. 2(t). We estimate that the SFE 
is 22 ± 8 % and that the break in the CMF occurs at a mass of 2.7 ± 1.3 Mq. Using the 
derived scaling factors for the IMF, we have also performed a Kolmogorov-Smirnov (KS) 
test between the scaled-up stellar IMF and the Pipe CMF. We find that the probability that 
the distributions are derived from the same parent population is 47%. 

For further comparison and to give a quantitative measure of how the inclusion of 
more information from the C^^O (1-0) emission effects the derived parameters, we have also 
performed the above analysis on the other three CMFs shown in Figure O Table [2] lists 
the derived parameters (mass scaling, SFE, and CMF break point) and their errors for each 
panel of the CMFs shown in Figure O Although Figures [5] (a) and (b) contain spurious and 
unrelated cores, we find that the derived parameters are similar to those determined for the 
other distributions. Considering each of the Pipe cores that have (1^0) emission as 

separate entities, as shown in Figure (c), we find that the CMF differs from the scaled IMF 
significantly for the highest masses. This is reflected in a low KS probability (7%) that the 
two are derived from the same parent population. Indeed, all three of these distributions have 
significantly lower probabilities (7-8%) of being derived from the same parent distribution 
as the stellar IMF compared to the adopted Pipe CMF (47%). 
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Discussion 



5.1. Derived core properties 

Because of the identification and selection methods used here, the individual core proper 
ties may differ slightly from those listed in our pr evious work (jAlves et al. 2007 : Muench et al. 



20071 : iLada et al.l l2008l : iRathborne et al.l l2008l ). However, the mean values for the radii, 
n(H2), and ant are identical to our previous work, i.e., mean radii of ~0.09pc, n(H2) of 
~7.3xl0'^ cm~'^, and ant of ~0.18 kms~^. Although the properties of the individual cores 
may differ, the ensemble properties remain unchanged: the cores appear to be mostly pres- 
sure confined entities whose properties are determined by the approximate balance between 
external and internal pressures coupled with self-gravity. Table [1] lists the derived properties 
for each core. 

By comparing the input ma sses to those derived for the cores when extracted using 
clfind2d, iKainulainen et al.l (120081 ) determined the overall uncertainty in the derived masses. 
There are three sources of error in the mass calculation. They are due to (1) the image 
noise, (2) the limiting extinction level considered by clumpfind, and (3) the degree of core 
crowding. The uncertainty in the derived masses due to the image noise can be deter- 
mined via the disp ersion in the ratio of the derived mass to true mass. The simulations of 
(120081 ) suggests that the derived masses have an uncertainty of ~ 25-30%. 



Kainulainen et al. 



Because the clfind2d algorithm only includes pixels down to a limiting extinction level (in 
this case an Ay of 1.2 ma gs), the core masses will tend to be systematically underestimated. 
Kainulainen et al.l (120081 ) show that the extent to which the derived masses are underesti- 
mated varies with core mass. For our clumpfind parameters, the simulations reveal that 
there is no correction necessary for cores with masses > 2 Mq. For cores with masses < 
2 Mq, however, the derived masses are typically 85% of the true mass. 

Masses may a l so be incorrectly determined if the cores are crowded and overlapping. 
Kainulainen et al.l (120081 ) defined a metric to describe core crowding which considers both 
the cores' relative separation and their diameter: f = mean core separation / mean core 
diameter. For cores within the Pipe Nebula, we find that their mean separation is ~ 0.38 
pc while their mean diameter is ~ 0.18 pc. The f-ratio is ~2.1, implying the crowding in the 
Pipe Nebula is minimal. Thus, the error in the cores masses will not be dominated by the 
effect of core crowding. 



- 15 - 



5.2. The Core Mass Function for the Pipe Nebula 



The derived CMF for the Pipe Nebula is shown as a probability density function in 
Figure [6] (solid curve). To generate the probability densi ty functions we use a Gaussian 
kernel of width 0.15 in units of log mass (ISilvermanI Il986l ) . Because this does not require 
discrete data bins, it more faithfully reproduces the detailed struc ture of the CMF c ompared 
to a binned mass function. Included on this plot is the CMF fromlAlves et al.l (120071 ) (dotted 



curve) and the stellar IMF for the Trapezium cluster (IMuench et al.ll2002l ) scaled in mass 
by a factor of ~ 4.5 (dashed line). The vertical dotted line marks the mass completeness 
limit of ~0.65 Mq. This limit gives the 90% mass completeness limit which reflects the core 
detections based on the extinction sensitivity, wavelet decomposition, and input clumpflnd 
parameters. The vertical dotted-dashed line marks the fldelity limit of ~1.1 M©, which 
corresponds to the mass above which the precise shape of the CMF is reliable. This limit 
is calculated by comparing the input and derived CMFs from the simulations and basically 
accounts for the fact that some cores, while reliably detected, have an incorrect mass due 
to the effects of core ove rlaps. Both of these limits were calculated from the simulations of 
Kainulainen et al.l (120081 ). 



Overall, th e shape of the CM F matches well with the scaled-up stellar IMF. Similar 
to the result of lAlves et al.l (120071 ). we flnd that the CMF is not characterized by a single 
power-law function. It can be described as a power-law shape above a mass of rrihreak ~ 
2.7 Mq, at which point there is a break from the single power-law and a clear flattening of 
the function towards lower masses. The presence of such a break above the completeness 
limit is signiflcant because it imparts a scale or characteristic mass to the CMF and enables 
a physically more meaningful comparison with the stellar IMF. The stellar IMF displays a 
similar break at mbreafc (IMF) ^ 0.6 Mq. The ratio of the two characteristic masses gives a 
direct measure of the SFE. As mentioned previously, we have performed a minimization 
between the stellar IMF and the Pipe CMF to give a quantitative measure of the SFE. 
Assuming that the dense cores in the Pipe will evolve to ultimately form stars, we flnd a 
characteristic SFE of ~ 22 ± 8% for the Pipe core population. 

While the mass completeness limit is ~ 0.65 Mq, all the cores included in the CMF 
below this mass are real cores and are associated with dense gas. Although the second break 
and apparent turnover of the CMF at M~ 0.4 Mq may be artiflcial due to the incompleteness 
of our sample, these features are tantalizingly close to a complementary second break and 
turnover in the scaled stellar IMF. The possibility that these features may be similar in both 
distributions potentially strengthens the connection between the CMF and the stellar IMF. 

Although the overall shapes of the Pipe CMF and stellar IMF are similar, the CMF 
does appear to fall off more steeply toward higher masses than the IMF. If this is signiflcant 
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it would suggest that the IMF is shghtly wider than the CMF and would require the SEE 
to be an increasing function of core mass in order to recover the stellar IME from the CME 
in a simple one-to-one transformation. However, given the increased uncertainty due to the 
sampling errors and small number statistics that characterize the higher mass bins, it is not 
possible to argue with any degree of confidence that the two distributions are different in 
this mass regime. 

The general similarity in the shapes of the CME and the IME holds significant impli- 
cations for understanding one of the key unsolved problems of star formation research, the 
origin of stellar mass. The fundamental nature of stars and their evolution has long been 
well explained by the theory of stellar structure and evolution developed in the last century. 
This theory has very successfully accounted for the mass- dependent luminosities, sizes and 
lifetimes of stars on the main sequence as well as their mass-dependent, post main-sequence 
evolution. However this theory is silent on the question of how stars get their masses in the 
first place. Because star formation research has established that stars form in dense molec- 
ular cloud cores, the similarity of the CME and stellar IME suggests that the IME derives 
its shape directly from the CME. If in fact the shapes of the CME and IME are the same, 
the derived SEE will apply to all core masses. This, in turn, would suggest that the IME 
results from a simple one-to-one transformation of cores into stars. Thus, it may be possible 
to trace the origin of the IME directly to the CME (modified by a star formation efficiency) 
and the origin of stellar masses directly to the origin dense core masses. 

However, the concept of a simple one-to-one transformation of cores into stars at a 
constant star formation efficiency is likely an oversimplification of the actual star formation 
process which is undoubtedly more complex. In particular, the more massive cores, that is, 
those cores in the cloud whose masses exceed the critical Bonnor-Ebert or Jeans mass are very 
likely to produce binary star systems. Indeed, the most massive cores may even fragment 
and produce small groups or clusters of stars. Thus, a strict one-to-one correspondence 
between cores and individual stars can not be preserved. However, the shape of the resulting 
IME could be retained, especially for higher stellar masses, if the SEE were to vary with 
mass, being higher for h igher mass cores, for example. Nonetheless, recent simulations by 



Swift fc Williamsl (120081 ) show that even when the internal fragmentation of cores in a CME 
is considered, the shape of the resulting IME is very similar to the shape of the input CME 
(apart from some details that may effect the very high and low mass ends in a small, but 
measurable way). Moreover, their results indicate that even if the SEE is not constant across 
the complete mass range, the resulting IMEs are not that different in shape from the original 
CME. Thus, while a one-to-one correspondence between cores and stars may not hold for all 
cores, the shape of the resulting IME is likely to be similar to the original shape of the CME 
and a characteristic or mean SEE can be measured by the ratio in the characteristic masses 
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of the two distributions. 



6. Summary and Conclusions 

In order to generate a more accurate CMF for the Pipe Nebula and thus enable a 
more detailed and meaningful comparison between the CMF and the stellar IMF, we have 
produced an improved census of the dense cores and their properties within the Pipe Nebula. 



Guided by the recent simulations of iKainulainen et al.l (120081 ) . we re-examined the ex- 



tinction map for the Pipe Nebula and extracted an improved and more reliable list of extinc- 
tion peaks and candidate dense cores. By systematically observing each of these extinction 
peaks in C^^O (1-0) emission, we have identified which peaks are associated with dense gas 
and which are associated with the Pipe molecular cloud. Moreover, we employed information 
about the kinematic state of the gas derived from our C^^O (1-0) survey to refine the iden- 
tifications of massive cores by distinguishing internal core structure from apparent structure 
caused by the overlap of physically discreet cores along similar lines-of-sight. The inclusion of 
spurious or unrelated cores, as well as the incorrect merging or separation of extinction fea- 
tures, can significantly affect the shape of the CMF constructed from observations. However, 
by combining measurements of visual extinction and molecular-line emission these effects can 
be minimized and more accurate measurements of core masses can be obtained. With more 
robust core identifications and masses we were able to construct an improved CMF for the 
Pipe Nebula. 

We confirm the earlier results which indicated a departure in the Pipe CMF from a 
single power-law form. We find the break point at a mass of ~ 2.7 ± 1.3 Mq, well above 
the mass completeness and fidelity limits for the observed core sample. Moreover, similar to 
Alves et al. (2007) we find that the overall shape of the CMF is generally similar to the stellar 
IMF except for a difference in mass scaling of about a factor of ~ 4.5. We interpret this 
difference in scaling to be a measure of the star formation efficiency (22 ± 8%) that will likely 
characterize the dense core population in the cloud at the end of the star formation process. 
These dense cores comprise an invaluable catalog for follow-up studies of the connection 
between core and star formation. 



We thank the referee for a thorough reading of the paper and the useful suggestions 
which have improved the paper. These observations were supported through NASA Origins 
grant NAG-13041 and the NASA Spitzer GO program (PID 20119) and supported by JPL 
contract 1279166. 
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Table 1. Core properties 
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Table 1 — Continued 
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*Core number from the catalog of lAlves et al.l (120071 ). Thirty three cores from the lAlves et al. 



( 120071 ) catalog are not included here because they either had no associated C^^O (1-0) emission or 
the C^^O (1-0) emission was not in the range 1.3 < vlsr< 6.4 kms~^. Because of the differences 
in the way the cores were extracted from the extinction image and the fact that we merged or 
separated t h e ext inction based on the associated C^^O (1-0) emission, many of the cores from 



Alves et al.l (120071 ) are listed multiple times in our catalo g. To be associated with a core in the 



current catalog, we require the cores of lAlves et al.l (120071 ) to have more than 20 pixels contained 
within the boundary of the new core. This criteria was necessary to prevent entries where the 
cores within the two catalogs overlapped only at the very edges. 

^Extinction peaks were merged within this core. 
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Table 2. Derived quantites from the minimization between the CMFs and IMFs in 

Figure [5l 



Panel 


Mass scaling 


SFE 


CMF break point 


KS probability 






(%) 


(Mo) 


(%) 


(a) 


3.6 ± 1.6 


28 ± 15 


2.2 ± 1.6 


8 


(b) 


3.8 ± 1.6 


26 ± 13 


2.3 ± 1.6 


7 


(c) 


4.5 ± 0.8 


23 ± 4 


2.7 ± 0.8 


7 


(d) 


4.5 ± 1.3 


22 ± 8 


2.7 ± 1.3 


47 
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Fig. 1. — Measured C^^O (1-0) central velocities (vlsr) toward the extinction peaks. The 
dotted vertical lines mark the range of velocities we define to determine which extinction 
peaks are associated with the Pipe Nebula. These values of 1.3 kms~^ and 6.4 kms~^, were 
determined from a Gaussian fit to the vlsr distribution. Of the 188 extinction peaks that 
are associated with C^^O (1-0) emission, we find that 158 of them have velocities in the 
range of 1.3 kms~^< vlsr< 6.4 kms~^. We consider these extinction peaks to be associated 
with the Pipe molecular cloud. All others are foreground and/or background molecular cores 
and are not included within the final list. 
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Fig. 2. — Three examples of how we determine whether adjacent extinction peaks are part 
of the same physical core or are unrelated. The left panels show the background-subtracted 
extinction images with contours as defined in § [2l The crosses mark the positions of each 
extinction peak identified. The right panels show the resulting cores after taking into account 
their velocity differences and Jeans lengths. Marked on these images are the central velocities 
determined from the C^^O (1-0) emission (V), the distance to the nearest extinction peak 
(D), and the Jeans length (Rj). In these images the color scale represents the area that is 
assigned to each core. In some cases we find that highly non-symmetric, complex extinction 
features have almost the same central velocity (top and middle panels). In other cases, 
however, we find that adjacent extinction peaks within the same large scale extinction feature 
can have very different velocities (lower panel). 
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Fig. 3. — The velocity difference (<5V) between each extinction peak and all others in a larger 
extinction feature plotted against the distance to the extinction peak's nearest neighbor (D). 
The filled circles mark the 6V and D for the extinction peaks that were merged into a single 
core. Forty one extinction peaks were found to be physically associated. These were merged 
into a total of 17 cores for inclusion in the final hst. 
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Fig. 4. — The 'cores- only', background-subtracted 2MASS visual extinction (Av) image 
( iLombardi et al.ll2006l ). Using the C^^O (1-0) emission toward each of the extinction peaks 
identified by clumpfind, we have determined which of the 201 extinction peaks are associated 
with dense gas, the Pipe Nebula, and if they are physically associated with any nearby 
extinction. Marked on this image are blue circles which correspond to the location and 
approximate extent of each of the 134 dense cores associated with the Pipe Nebula. The green 
plus symbols mark the extinction peaks that are associated with foreground/background 
molecular clouds. The red crosses mark the extinction peaks that had no detectable C^^^O (1- 
0) emission. 
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Fig. 5. — Derived CMFs shown as binned histograms. The four panels correspond to dif- 
ferent core samples when taking into consideration the C^^O (1-0) molecular line emission. 
Note that the shape of the mass function changes considerably when we use the C^^O (1-0) 
emission to guide the core extraction from the extinction image (i .e. between panels (a ) and 
(d)). The dashed line corresponds to the scaled field star IMF of iMuench et al.l (120021 ). For 
each panel we determine, via a minimization, the offsets between the CMF and the IMF. 
To accurately determine these scaling factors, we minimize the between the distributions 
only for masses above the completeness limit. The derived parameters are summarized in 



Table [2l The vertical dotted line marks the mass completeness limit (IKainulainen et al. 



20081 ). Included on the histograms are the errors for each bin (calculated as the square root 
of the number per bin). 
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Fig. 6. — The improved CMF f or the Pipe Nebula as a probabihty density function (solid 



curve) compared to the CMF of lAlves et al.l (l2007f) ( dotted curve). The dashed line corre- 
sponds to the field star IMF of iMuench et al.l (120021 ) scaled up by a factor of ~ 4.5. We 



interpret this difference in scaling to be a measure of the star formation efficiency (22 ± 
8%) that will likely characterize the dense core population in the cloud at the end of the 
star formation process. The vertical dotted line mar ks the mass completene ss limit, while 
the vertical dot-dashed line marks the fidelity limit (IKainulainen et al.ll2008l ). We confirm 
earlier results that suggested that the Pipe CMF departs from a single power-law with a 
break at ~ 2.7 ± 1.3 Mq. 



